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Since the beginning of the efforts to produce and understand quark matter large changes 
developed in the ideas of description of this matter. In the present paper we summarize 
some aspects of this development. 

1. Evolution of ideas on quark matter 

Since the beginning of the quest of quark matter deep changes have developed in our 
concept in understanding its properties. In Table [T] the ideas are collected, which domi¬ 
nated at the beginning of quark matter research. The present status of these ideas are also 
displayed. Although theoretical models supported the early concept, but the experimental 
data forced to change these speculations step by step. Thus the picture of a static, weakly 
interacting gaseous quark-gluon plasma was substituted by the more dynamic, strongly 
interacting fluid-like quark-gluon matter, where the ingredients gained an effective mass 
during the intense microscopical interactions. 


Early times 

Recent status 

very high quark density; 

not too high quark density; 

vanishing coupling constant; 

large effective coupling constant; 

zero mass for quarks and gluons; 

hnite effective mass for quarks and gluons; 

weakly interacting constituents; 

strongly interacting constituents; 

static, homogeneous gas; 

dynamic, liquid-like system 

very high temperature; 

temperature near to critical; 

vanishing sound velocity; 

hnite sound velocity; 

long lifetime; 

short lifetime; 


Table 1 

Concept of Quark Matter in the early times [^J and recently, considering the microstruc¬ 
ture (upper part of the table) and the bulk properties (lower part of the table). 
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2. Time evolution of a heavy ion reaction 

As soon as the static picture of quark gluon plasma was substituted by a dynamic 
scenario with short lifetime, the different pre-QGP and post-QGP stages became as im¬ 
portant as the QGP state itself. As the dominance of a static plasma state has been lost, 
the dynamical steps and the intermediate phases (’matters’) became the target of the 
theoretical investigations. Table 2 displays the emerged scenarios appearing in heavy ion 
collisions. 


type of event 

microscopic dynamics 


color glass condensate (?) 

gluon saturation [E] 


evolution of plasma 

expands and cools. 



undergoes through a set of phases [E] 


prehadronization stage 

quarks and antiquarks gather effective mass 

[m 

hadronization 

quarks and antiquarks coalesce into hadrons 

[Q 


Table 2 

Time evolution of a heavy ion reaction. 


In a heavy ion collision not a single sort, well defined type of matter is formed, but a 
set of different types of matter is created in the consecutive stages of the reaction. 

3. Evolution of ideas on hadronization 

In case of a quasi-static quark gluon plasma with long lifetime, there should be enough 
time for slow hadronization, which is driven by energy and momentum transformation at 
the macroscopical level. This leads to the application of thermodynamics. 15 years ago 
we were very much confident to expect a thermodynamical first order phase transition. As 
the time scale of the reaction became shorter, the idea of mixed phase transition appeared, 
handled with thermodynamical methods. However, as it turned out that the evolution 
of a heavy ion reaction is much faster, the hadronization ideas too had to be adjusted to 
this faster scenario. Finally, so fast hadronization processes must have been considered, 
that the macroscopical description was substituted by microscopical one and the slow 
thermodynamical transformation has been substituted by the fast process of coalescence 
of constituent quarks. 

4. The coalescence models 

The concept of quark coalescence has been evolved through different steps. At first the 
idea of quark number conservation appeared and became a very useful tool [Ej. Then 
the linear vs. nonlinear coalescence descriptions were investigated [ Later on the 
description was improved and the direct momentum conservation at the microscopical 
level was considered in the coalescence processes [EIEIE]. Finally mass conservation has 
also been studied to produce light hadrons from heavy constituent quarks [CHI- 
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4.1. Coalescence hadronization with constituent quark number conservation 

For the description of hadrons the constituent quark model was quite successful. In this 
approach the quarks are dressed, with mass of approximatively 300 MeV. The question 
arises, that at what stage of hadronization is this effective mass created? In the ALCOR 
model [Ej it is assumed that this effective mass is created in the very last stage, in the 
prehadronization stage of the evolution of fireball. 

The main assumptions of the ALCOR model: 

A) At the beginning of hadronization the quarks are dressed constituent quarks. 

B) These quarks coalesce to form the hadrons. 

C) The number of different quarks and antiquarks is conserved during hadronization. 

D) The effective mass of gluons is much higher than that of quarks near the critical 

temperature. Thus the gluon degree of freedom is neglected in this late period of 
the heavy ion reaction. 

E) The number of a given type of hadrons is proportional to the product of the numbers 

of different quarks from which the hadron consists: 

As,(pfc) = Csxijk) ■ {h ■ Ni) ■ {bj ■ Nj) ■ • Nk) (1) 

where NB,(ijk) is the number of produced baryons from quark i, j, fc, and the equation 
for quark number conservation in the hadronization ]Sff^^^d.^orimatter _ _ jy. 

determines the normalization coefficients. 

With these assumptions one gets the relative numbers of particles in good agreement 
with experimental data, see Table 01 



ALCOR model 

STAR data 

Ref. 

h- 

280 

280 ± 20 

[EH 

K-/71- 

0.159 

0.161 ±0.002 

m 

K+/K- 

1.091 

1.092 ±0.023 

m 

■p- /p+ 

0.66 

0.71 ±0.05 

[El 

A/A 

0.72 

0.71 ±0.01 

m 


0.80 

0.83 ±0.04 

m 

E-/h- 

0.0091 

0.0077 ±0.001 

m 

^/K*° 

0.42 

0.49 ±0.12 

[ESI 


Table 3 

Hadron production in Au+Au collision at ^/s = 130 AGeV from the ALCOR model and 
experimental data from the STAR Collaboration [ ITTl IT^ OIT^ ITHl ITT)] . 
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4.2. The simple quark counting 

Realizing [El that the coalescence probabilities for a hadron and its antiparticle are the 
same: Ch = C-j ^, very transparent relations can be obtained for the ratios of the different 
mnltiplicities: 


A 

A 






( 2 ) 


These predictions agree again with the experimental data, see Table Eland Figure H 


Ratios 

STAR 

SPS 

A+/A- 

1.092 ±0.023 

1.76 ±0.06 

{A/A}/{p/p} 

0.98 ±0.09 

2.07 ±0.21 

{WS}/{A/A} 

1.17±0.11 

1.78 ±0.15 

{n/G}/{S/S} 

1.14 ±0.21 

1.42 ±0.22 


Table 4 

Compilation of experimental particle ratios obtained at RHIC energy y/s = 130 AGeV [ 
EH] and SPS energy 17.3 AGeV [ IT71 ITHl ITUl EH] . 


2.25 
2 

1.75 
1,5 

1.25 
1 

0,75 

0.5 

0.25 

0 


1 Pb + Pb at SPS 


1 Au+Au at RHIC 

f. 


+ 

i. 


h! a/a 

K' p/p A/a 


n/Q 


Figure 1. Graphical display of the particle ratios at mid-rapidity in Au-|-Au collision at 
RHIG [EH] and in Pb-|-Pb collisions at SPS [EH EH EH 1201 • 
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4.3. Charge fluctuation in a quark-antiqnark system 

The charge fluctuation [ini is characterized by the ratio 


D 


^ <SQ^> 

< > ’ 


( 3 ) 


where < dQ‘^ > is the average of the charge fluctuation and < Nch > is the average value 
of the charge multiplicity. The value of D is D ~ 3 for hadron gas in equilibrium, and 
D ^ 1 for the quark gluon plasma, The measured D value is near to 3. However, in the 
constituent quark coalescence scenario (ALCOR) one can also expect a D ~ 3 value. 


4.4. Coalescence with constituent quark momentum conservation 

The idea of quark coalescence can be applied at relatively large momenta and the hadron 
spectra at high-pj^ can be investigated successfully [ El |Hl E) • In this kinematic region, 
especially at > 2 GeV, the mass of the light and strange quarks can be neglected, 
Pquark » ^quarki and the Calculations are simplified. Assuming that the most dominant 
kinematic window is overwhelming, then 


Pa=Pb = Ph/2 = P/2 


( 4 ) 


The expression for the meson emission spectrum [E]: 


.dAT, 


M 


(PP 


= da 


P ■ u{r) 
(27r)3 




P. 


fair; -) fb{r-, -). 


E 


A similar expression is valid for the baryon emission spectrum: 

JN, 


B 


d^P 


= da 


P ■ u{r) 
(27r)3 


P. 


P. 


fair; —) fb{r-, —) fc{r-, —). 


( 5 ) 

( 6 ) 


These two equations lead to an interesting consequence: since the momentum distribu¬ 
tion of constituent quarks is expected to behave as an exponentially decreasing function 
of transverse momentum, there are more partons at P/3 than at P/2. 

fqiP j) > fq{r-, (7) 

From that follows, that we have more baryons than mesons in the momentum range 
where hadrons are produced by coalescence. This explains the surprising momentum 
dependence of the observed proton to meson ratio [EIEIEI. This description also implies 
the conservation of quark and antiquark numbers in the hadronization process. 


4.5. The effective mass of quarks 

Mass distribution of quarks 

Since the effective mass of quarks are determined by their interaction with the neigh¬ 
borhood, and the matter in the fireball is not smooth and homogeneous, the granulated 
structure of the fireball is reflected in the effective mass distribution of the quarks. 

The consequences of a possible quark mass distribution 

_ _M_ / /J. I ttT 

p{m) = Ne TcV 

with p = 0.25 GeV and Tc = 0.26 GeV, are discussed in the next paragraph [nni- 


( 8 ) 
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Figure 2. The product of the two quark 
mass distributions p{m) ■ p{mT^ — m) 
in eq. m- The maximum is at m = 
m^/2. 


Figure 3. Product of three mass distributions 
in eq. (USD- The horizontal axes are mi and 
m 2 , the mass m 3 is dehned as m 3 = tub — 
(mi + m2). 


4.6. The coalescence expression for distributed quark mass 

Using the covariant coalescence model of Dover et al, [123, the spectra of hadrons 
formed from the coalescence of quark and antiquarks can be written as 



(9) 


Here Fh{ph]Xh) is an eight dimensional distribution (Wigner function) of the formed 
hadron, and da denotes an infinitesimal element of the space-like hypersurface of hadron 
production. 

Assuming that meson production is homogeneous within the reaction volume, the source 
function for the a + b ^ M meson production, Fm{pm), is the following: 



( 10 ) 



( 11 ) 


The parameters om and reflect properties of the hadronic wave function in the mo¬ 
mentum representation convoluted with the formation matrix element. 

Assume that Pc is so small, that partons with practically zero relative momentum form 
a meson. 


Pa = Pb = Pm/2 ma + mb = tum 


ma + mb = ruM 


( 12 ) 


This leads to the meson coalescence function 


Cm = c<m- i(Pa - Pm/ 2) ■ (llPb - Pm/ 2) ■ i5(m„ + m, - TOm) 


(13) 


Thus we arrive at the following meson distribution function: 



p(ma) ■ fq(pt/2;ma) ■ p(mb) ■ fb{pt/2,]mb) . 


( 14 ) 
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Figure 4. Measured [ |221 and calculated 
7r° and 7r“ transvers momentum spectra 
in central Au+Au collisions at a/s = 200 
AGeV. 


Figure 5. Measured [ |221 and calcu¬ 
lated transverse momentum spectrum for 
antiproton production in central Au-|-Au 
collisions at y/s = 200 AGeV. 


Similar argumentation leads to the three-fold coalescence expression. In the place of 
eq. m we get: 

Pi = P2 = Ps = Pb/3 mi + m2 + ms = niB (15) 

and the baryon source function becomes 

rTTLpr I'Tflpr rTTlpr 

FsiPt) = CiB dmi / dm 2 / dm 3 S{mB - {mi + m 2 + m 3 ) ■ 

JO Jo Jo 

■p{mi) /g(pi/3,mi) ■ p{m 2 ) fq{pt/3,m2) - pirns) fqipt/S^ms) . (16) 

The maximum contribution to the coalescence integral is obtained from the equal mass 
part of the mass distributions. 

The calculated meson and baryon source functions and their space-time integrals in 
eq. © yield theoretical spectra comparable to experimental data [EB]- In Figure IH we 
display our results for pion production together with the experimental data on 7 r° and 
vr"'' production at mid-rapidity in Au-|-Au collisions at y/s = 200 AGeV. Figure El shows 
calculation and data for antiproton. The agreement in the 0.5 < pt < 5.0 GeV momentum 
window is astonishing and proves the validity of the coalescence scenario. 

Besides the hadron production on a logarithmic scale, we display the p~/ ti~ ratio in 
linear scale on Figure 6 . The theoretical result is able to follow the tendencies of the 
anomalous antiproton to pion ratio. 
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Figure 6. p/n ratio as a function of pt in central Au+Au collisions at = 200 AGeV. 
Experimental data points are from Ref. [I20j. 


4.7. The v2 paradox 

The meson and baryon elliptic flow at high pt show a peculiar behaviour, which is 
understandable in the framework of parton coalescence dynamics [Ell- 

Assuming the following form for the asymmetric quark momentum distribution: 


dNq 

pj_dp±d^ 


1 dNq 
2tt p±dp± 


(1 + 2v2,q{q±)cos{2^)) 


(17) 


one gets in the coalescence dynamics for U 2 -C 1 the hadron elliptic flow: 


n2,M(p±) ~ ‘2*V2,q{p±/2) (18) 

V2,b{p±) ~ 3*U2,q(p±/3) 

This behaviour of the V 2 ,m and U 2 ,b parameters predicted by the coalescence hadroniza- 
tion dynamics are clearly observable in the experimental data (see Fig. El). 

The behaviour of flow asymmetry parameter 


V2,u — V2,d — V2,s 

implies that the collective flow for quarks of all flavour is the same. 


( 19 ) 
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Figure 7. Measured baryon and meson elliptic flow (flgure is from Ref. [ESI)- 


5. Conclusion 

A large number of experimental facts is in agreement with the assumption that in the 
prehadronization stage 

• the quark matter consists of constituent quarks and antiquarks with effective mass; 

• the collective flow for quarks of all flavor is the same; 

• the hadronization proceeds via coalescence mechanism; 

• the numbers of quarks and antiquarks are conserved during hadronization. 

Further 

• the initial stage of heavy ion collision is strong color held (gluon) dominated; 

• the final stage of heavy ion collision is quark dominated. 
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